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Abstract Tracking the sign of fluctuations governed by the (1 + 1)-dimensional Kardar-Parisi-Zhang (KPZ) 
universality class, we show, both experimentally and numerically, that its evolution has an unexpected link to 
a simple stochastic model called the renewal process, studied in the context of aging and ergodicity break¬ 
ing. Although KPZ and the renewal process are fundamentally different in many aspects, we And remarkable 
agreement in some of the time correlation properties, such as the recurrence time distributions and the persis¬ 
tence probability, while the two systems can be different in other properties. Moreover, we And inequivalence 
between long-time and ensemble averages in the fraction of time occupied by a specific sign of the KPZ- 
class fluctuations. The distribution of its long-time average converges to nontrivial broad functions, which are 
found to differ significantly from that of the renewal process, but instead be characteristic of KPZ. Thus, we 
obtain a new type of ergodicity breaking for such systems with many-body interactions. Our analysis also 
detects qualitative differences in time-correlation properties of circular and flat KPZ-class interfaces, which 
were suggested from previous experiments and simulations but still remain theoretically unexplained. 
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Keywords Growth phenomenon ■ Scaling laws • KPZ universality class • Renewal theory • Stochastic 
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1 Introduction 

The Kardar-Parisi-Zhang (KPZ) universality class ll^ l2ll27l[8l is a prominent nonequilibrium class, ruling di¬ 
verse kinds of nonlinear fluctuations in growing interfaces Il26ll2ll46l . driven particle systems Il27l [8l. fluctuat¬ 
ing hydrodynamics 1411 . and so on. Particularly noteworthy are recent analytical developments on the (1 -f 1)- 
dimensional KPZ class, which have exactly determined a number of its statistical properties on the solid math¬ 
ematical basis ll27]l8]l . Specifically, for (1 -h 1)-dimensional KPZ-class interfaces, the interface height h{x,t), 
measured along the growth direction at lateral position x and time t, grows as 

/r(x,f) + (1) 

with parameters Voc and F, a rescaled random variable X, and j3 = 1 /3 being the characteristic growth exponent 
of the (1 + 1)-dimensional KPZ class Il26l l2). Then the recent analytical studies ll27l [8l have consistently 
shown that z(-^T) exhibits one of the few universal distribution functions, selected by the choice of the initial 
condition, or equivalently the global shape of the interfaces. For example, circular interfaces grown from 
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a point nucleus show the largest-eigenvalue distribution of random matrices in Gaussian unitary ensemble, 
called the GUE Tracy-Widom distribution, while flat interfaces on a linear substrate show the equivalent for 
Gaussian orthogonal ensemble. This implies that the KPZ class splits into a few universality subclasses. They 
are also characterized by different spatial correlation functions, whose exact forms are also known analytically 
iT/l ISl. These results for the circular and flat subclasses were corroborated by direct experimental verifications 
using growing interfaces of turbulent liquid crystal (LC) Il47ll49ll48l . 

In contrast to these clear characterizations of the distribution and the spatial correlation, analytical results 
on the temporal correlation remain limited, hence challenging II12II24I[T31[T^ . The LC experiment ll48l and 
numerical simulations 1125 114011451 1711 showed that the temporal correlation is also different between the circular 
and flat subclasses. Firstly, the two-time correlation function 0 

C(f,fo) = {h{x,t)h{x,to)) - (/r(x,f))(/r(x,fo)) 

(r2fof)'/3c,e,(tAo) (2) 

was shown to decay, in its rescaled form, as Cres(f/G) (t/to) ^ with A = 1 for the flat case 14811251 171 
and A = 1 /3 for the circular case Il48ll40ll45l l7l. The latter implies lim,_>„c C(f, g) > 0, i.e., correlation remains 
strictly positive, forever, in the circular case. Secondly, the persistence probability P± (f, to) was also measured, 
which is defined here by the probability that the fluctuation 5h{x,t) = h{x,t) — {h{x,t)) at a fixed position x 
never changes its sign (from the one denoted by the subscript) in the time interval This quantity was 

found to show a power-law decay 

■P±(fTo)(3) 

with exponents w 1.35 < w 1.85 for the flat interfaces II25II48PI and w 0.8 for the circular 

ones Il40ll48ll45l . The latter implies divergence of the mean persistence time /“P±(f,fo)dt in the circular case. 

It has been shown that such a divergent mean leads to anomalous dynamics such as non-ergodicity, anoma¬ 
lous diffusion, aging, and population splitting II18II33II38II39I . Therefore, the above observations on the circular 
KPZ subclass imply that it may also be understood in this line of research. Ergodicity is a basic concept in 
statistical physics and dynamical systems. It guarantees that time-averaged observables obtained by single 
trajectories converge to a constant (ensemble average) as time goes on. Ergodicity breaks down when, e.g., 
the phase space consists of mutually inaccessible regions, because then single trajectories are unable to cover 
the whole phase space. However, in 1992, Bouchaud ||3 proposed another situation of ergodicity breaking, 
where the phase space is not split, but trajectories undergo long and random trapping. If the mean trapping 
time diverges, trajectories cannot sufficiently explore the phase space, however long they do. Interestingly, 
in this situation named weak ergodicity breaking (WEB) ||3, certain time-averaged observables such as the 
time-averaged diffusion coefficient Il20ll35l fn [3^l34l and the fraction of time occupied by a given state OOl 
fT8ll32ll33l do not converge to their ensemble average, but themselves become well-defined random variables, 
described by characteristic distribution functions in simple cases. The existence of such an asymptotic broad 
distribution for time-averaged quantities is usually regarded as a defining property of WEB. Experimentally, 
single-particle observations have indeed shown relevance of WEB in macromolecule diffusion in biologi¬ 
cal systems 05 1 1123 Il50ll44ll3 1 11341 and in blinking quantum dots Il29l l5l l43ll42l . However, it remains unclear, 
both theoretically and experimentally, how useful these developments on WEB are to characterize many-body 
problems such as the KPZ class. Therefore, it is a challenging and important issue to clarify if WEB occurs 
in KPZ, and if yes, characterize the WEB of the KPZ class. 

* Throughout the paper, (• • ■) denotes the ensemble average defined over infinitely many realizations. It is independent of 
X because of the translational symmetry. Therefore, for evaluation, we can take averages over positions too to achieve better 
statistical accuracy, without changing its mathematical definition. 

^ In the literature, our definition of P±(t,to) based on the ensemble average is sometimes called the “survival probability”, 
in which case the term “persistence probability” is reserved for the probability that the sign of h(x,t) — h{x,to) is unchanged 
l4l . However, in the present paper, we define our persistence probability by the ensemble average (unless otherwise stipulated), 
following earlier studies of direct relevance. 

^ holds when the underlying KPZ nonlinearity ^ (V/t)^ is positive, otherwise the order is reversed 1^ . Note also 

that the estimates obtained in (25] are somewhat different from those from the LC experiment 1481 . We believe this is because 
of the different choice of the reference time to: while it was taken to be right after the initial condition in 1251 . times in the 
asymptotic KPZ regime were used in the LC experiment (48). This is also the choice of the present paper; therefore, here we 

refer to the estimates from the LC experiment, 0^^ ~ 1-35 and ps 1.85, as the values of the persistence exponents for the 
flat KPZ subclass. 
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Fig. 1 Waiting-time distributions (ccdfs) /7±(T;r) against zjt for the circular interfaces, obtained at different t in the LC exper¬ 
iment (main panels) and the Eden model (insets). The dashed lines are guides to the eyes indicating exponent —0.8. The same 
set of colors/symbols and t is used in both panels. The ordinates are arbitrarily shifted. 


2 Our approach 

To investigate possible relationship between WEB and KPZ, we construct a dichotomous process from height 
fluctuations of KPZ-class interfaces, a{x,t} = sign[5/t(x,t)], which is then regarded as a time series. Such 
dichotomization has recently been used to characterize time-correlation properties of interactions on lipid 
membranes 15211531 and of turbulence I22II21I . successfully. The constructed process is compared with a 
theoretically defined dichotomous process, arguably the simplest and best-studied one, namely the renewal 
process (RP) lfT5l l9lfT8l. RP consists of a single two-state variable, which switches from one to the other state 
after random, uncorrelated waiting times generated by a power-law distribution: 

/ T \ 

p { t ) = Prob[waiting time > t] = I — , (t > Tq). (4) 

V^o/ 


This model shows WEB and aging for 0 < 0 < 1 lfT5l l9l [T8ll39l . 

Concerning KPZ-class interfaces, we use the experimental data of the circular and flat interfaces ob¬ 
tained in Refs. I49II48I (LC turbulence): for the circular (or flat) case, total observation time was Ttot = 
30.5 s (63 s), time resolution was Tres = 0.5 s (0.35 s), and N = 955 (1128) realizations were used, respec¬ 
tively. We also analyze newly obtained numerical data for circular interfaces of the off-lattice Eden model 
||45] (Ttot = 5000,7(es ~ I,N = 5000) and flat interfaces of the discrete polynuclear growth (dPNG) model 
(Ttot = 10^, Tres —0.1,N— 10"*). Eurther descriptions of the systems and parameters are given in [Appendix A| 


3 Results 

3.1 Circular interfaces 

Eirst of all, we stress that RP is far too simple to fully describe KPZ, because RP is a two-state model without 
even spatial degrees of freedom and has uncorrelated waiting times. We nonetheless measure the waiting times 
between two sign changes of 5h, first for the circular interfaces, for which we anticipate relation to WEB as 
discussed above. More specifically, we define the waiting-time distribution p±{z\t) by the probability that the 
sign renewed at time t (changed to the subscripted one) lasts over time length t or longer, hence p± (t; f ) is the 
complementary cumulative distribution function (ccdf). Eigure[T] shows the results for both the LC experiment 
(main panels) and the Eden model (insets). Remarkably, in both cases we And a clear power law as described 
in Eq. (|4]i with exponent 6 = 0.8, while the cutoff Tq is out of the range of our resolution. 

This similarity to RP leads us to compare further statistical properties between the two systems. Eirst we 
focus on the forward recurrence time Tf(f ), defined as the interval between time t and the next sign change, as 
well as the backward recurrence time Tb(t), which is the backward interval from t to the previous sign change 
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(a) circular, forward, + sign (b) circular, forward, - sign 





Fig. 2 Distributions (pdfs) of the forward (a,b) and backward (c,d) recurrence times, Zf and Tb. respectively, for the circular 
interfaces at different t in the LC experiment (main panels) and the Eden model (insets). Tf and Tb are rescaled by t. The black 
lines indicate RP’s exact results, Eqs. 0 and 0, with 6 = 0.8. The data are normalized so that they have the same statistical 
weight as RP’s exact results in the range covered by their abscissa. The same set of colors/symbols and t is used in all panels. 
Note that data at f = 25 s are not shown for Tf because the remaining time is then too short to measure the distribution of Tf. 


lll8lll4lfT5l . For RP, Dynkin 11411181 derived exact forms of the probability density function (pdf) of Tf(f ) and 
Tb(f) as follows, for 0 < 0 < 1: 


pdf(ff) = 


sin n9 1 

n ff®(l + ff)’ 


pdf(fb) = j3i_0,e(fb) = ®(l-fb)® 

n 


(5) 

( 6 ) 


with Tf = Tf/f, fb = Tb/f and I5a,b{x) denoting the pdf of the beta distribution. Although their derivation essen¬ 
tially relies on the independence of waiting times in RP, a feature not shared with KPZ, we hnd, as shown in 
Fig. 121 that both experimental and numerical results for the circular interfaces precisely follow RP’s exact re¬ 
sults indicated by the solid lines (except finite-time corrections). Note that the persistence probability P± {t, to) 
considered in Eq. (O actually amounts to the ccdf of Tf(fo), i.e., P±{t,to) — pdf('^f(fo))dTf. This indicates 
that the functional form of the persistence probability, which is usually intractable for such spatially-extended 
nonlinear systems lU, seems to be given by RP’s exact result (O in the case of the circular KPZ subclass. We 
also remark that the explicit dependence of the pdfs on t indicates the aging of the system. 

In contrast to this agreement, we also find statistical properties that are clearly different between the two 
systems. The occupation time T+, i.e., the length of time spent by the positive sign, is a quantity well-studied 
in two-state stochastic processes OOII18I and in more general scale-invariant phenomena (see, e.g., Illll l. It is 
simply related to the time-averaged sign a = (1/T) Jq a(x,t)dt by a = 2T+ jT — 1. For RP with 0 < 0 < 1, 
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Fig. 3 Time-averaged sign distribution (a) and sign coiTelation function (b) for the circular interfaces, (a) Pdf of the time- 
averaged sign O' = (1/T)fg cj(x,?)df for the LC experiment (T = 7.5s, 15s,30.5s for circles, squares, and diamonds, respec¬ 
tively) and the Eden model (T = 5000), compared to the Lamperti distribution 0 with 6 = 0.8 (dashed line). The gray vertical 
line indicates the ensemble-averaged value (a) = —0.021. The existence of the broad asymptotic distribution is a direct evidence 
of WEB in the circular KPZ subclass, (b) Correlation function of sign, Csign(t,to) = (o’(x,t)cF(x,to)), measured at different to for 
the LC experiment (main panel) and the Eden model (inset). The dashed and dotted lines indicate exponents —A = —1/3 and 
— 6 = —0.8, respectively. 


Lamperti ll^ showed that it does not converge to the ensemble average, but remains stochastic even for 
L —oo, with its pdf derived exactly as follows O0II18II : 


pdf(a) 


(2sin7re/7r)(l — ^ 

(l+a)2® + (l-a)2® +2(cos7r0)(l—a^)®' 


(7) 


This distributional behavior of the time-averaged sign is a clear evidence of WEB in RP. The corresponding 
pdfs obtained at different T for the circular KPZ interfaces [Fig. |2a) symbols] indeed indicate an asymptotic 
broad distribution, demonstrating that a remains stochastic and does not converge to the ensemble average 
(a) = —0.021 (determined by the GUE Tracy-Widom distribution) shown by the gray vertial line in Fig.[3[a). 
This demonstrates that KPZ indeed exhibits WEB, at least for the circular case. On the other hand, the found 
distribution is clearly different from the Lamperti’s one for RP with 6 = 0.8 (black dashed line). We find 
instead a nontrivial distribution universal within the circular KPZ subclass, as supported by good agreement 
between experiments and simulations (symbols and turquoise solid line). 

Another quantity of interest is the correlation function of sign, Csign(f,fo) = (<7(x,f)a(x,fo)). This can be 
expanded by the generalized persistence probability P±,n(Lfo )5 i-C-^ the probability that the sign changes n 
times between Iq and t (hence P± o{t,to) =p±{t,to)y- 


Csign(Lfo) = ^7^(0) ^ (-l)"-P±,H(Lfo), 
± «=0 


( 8 ) 


where P± (fp) denotes the probability that fluctuations at fp take the sign indicated by the subscript. For RP with 
0 < 0 < 1, one can explicitly calculate the infinite sumof Eq. ® and obtain Csign(f,fo) — L±7’±(fo)7’±,p(f,fo) 

(f/fp)“® ITSll . In contrast, for the circular KPZ interfaces, we And that Csign(f,fo) decays as Csign(f,fo) ~ 
with 1 = 1/3 [Fig.[3b)] (see also footnote |4]i, in the same way as the rescaled correlation function Cres(t/fo) 
does [Eq. (|2]l]. Since the relation ® holds generally and P±p(f,fp) = pdf(Tf(fp))dTf is alike, the differ¬ 
ence from RP should stem from P±.„>i (f,fp), which encodes correlation between waiting times. For RP, one 
can show 


P±,n>\{t0 + 


Af-2®+i for At < to, 
At~^ for At ^ to. 


(9) 


with to, At ^ Tp (see [Appendix B[ ). Now, for the circular KPZ interfaces, the results in Fig. |4] show that the 
long-time behavior seems to be consistent with that of RfQ, but the short-time behavior for odd n shows faster 
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For some quantities the asymptotic decay is only reached hy the numerical data, obtained with longer time. 
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(a) LC circ. at = 4s, + sign (b) LC circ. at = 4s, - sign 






Ml tQ 


Ml tQ 


Fig. 4 Generalized persistence probabilities + At,to) of the circular interfaces, measured for the positive (a,c) and 

negative (b,d) fluctuations (sign at tg is used) in the LC experiment (a,b) and the Eden model (c,d). The dashed and dotted lines 
indicate exponents —29 + 1 = —0.6 and —9 = —0.8, respectively. The same set of colors/symbols and tt is used in all panels. 


decay than that of RP (after the initial growth, which occurs at Zif < n*/® To for RP; see [Appendix B[ l. In other 
words, P±^„(to + At,to) has heavier weight in the short-time regime for odd n. This difference from RP gives 
nontrivial contribution to the sum in Eq. which is absent for RP. We consider that this is how the different 
behavior of the correlation function arises, which captures, for KPZ, the characteristic time correlation of the 
(non-binarized) KPZ-class fluctuations. 


3.2 Flat interfaces 

Now we turn our attention to the flat interfaces. Figure |5] shows the waiting-time distribution (ccdf) p±{z\t) 
for the flat LC experiment (main panels) and the dPNG model (insets). At short waiting times, we identify 
power-law decay with exponent —0.8. This exponent seems to be different from —2/3 previously observed 
for a related quantity in the KPZ stationary state 125 iR . but is identical to the one we found for the circular 
interfaces (Fig. [T]l. For the flat interfaces (Fig.|5]l, however, this power law is followed by another one with 
larger (in magnitude) exponent for longer waiting times, which now takes different values between positive 
and negative fluctuations. The measured exponents do not seem to reach their asymptotic values within our 
observation time, increasing gradually with T, but they are clearly asymmetric with respect to the sign, in sharp 
contrast with the exponent for the shorter waiting times or for the circular interfaces. Moreover, p± (t; f ) taken 
at different f overlaps when it is plotted against z/t (Fig.|5]l. This indicates that the value of T separating the 

^ The exponent measured in 1^ was about the persistence probability of the sign of Ah{x,t) — Ah(x,tg) with Ah(x,t) = 
h{x,t) — fg h(x,t)dx/L, whereas we measure here the waiting-time distribution, or the persistence probability of the sign of 
Sh{x,t) = h(x,t) — {h{x,t)) with the condition 8h(x,tg) = 0. Although both probabilities concern the first return to zero, the 
different definitions may lead to different exponent values. 
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Fig. 5 Waiting-time distributions (ccdfs) p±{T;t) against f/t for the flat interfaces, obtained at different t in the LC experiment 
(main panels) and the dPNG model (insets). The dashed and dotted lines are guides to the eyes, indicating exponents labeled 
alongside, though for the dotted lines we expect larger asymptotic exponents (see text). The same set of colors/symbols and t is 
used in both panels. The ordinates are arbitrarily shifted. 



(c) flat, backward, + sign (d) flat, backward, - sign 



Fig. 6 Distributions (pdfs) of the forward (a,b) and backward (c,d) recurrence times, Tf and Tb, respectively, for the flat interfaces 
at different t in the LC experiment (main panels) and the dPNG model (insets). Tf and Tb are rescaled by t. The black lines indicate 
numerical results for the 2-step RP at ? = 10^. The data are normalized so that they have the same statistical weight in the range 
covered by their abscissa. The same colors/symbols correspond to the same t in all panels. 


two power-law regimes is not constant, but grows with t, showing the aging property of the waiting-time 
distribution. 















































These results on the flat-KPZ waiting-time distribution lead us to introduce a variant of RP with two 
power-law regimes, called hereafter the 2-step RP model0: 


p±(T;f) = < 




for To < T < f, 
for f < T. 


( 10 ) 


We then solved it numerically with 0 = 0.8, 0+ = 1-2, 0'_ = 1 -5 (as observed experimentally) in the following 
way: First, the initial sign was chosen to be either -|- or — with the equal probability. The first waiting time 
was generated according to p±{x\t) — (t/to)~®± (t > To), because Eq. (fTOl) is invalid for t — 0. Subsequent 
waiting times were generated by Eq. (1101) . until the time (cumulative sum of waiting times) exceeds the 
recording time Ti^t- We sampled 10^ independent realizations to investigate statistical properties of this 2-step 
RP model. 

Now we compare this 2-step RP model and the flat KPZ interfaces. Pigure|6]shows the forward and back¬ 
ward recurrence-time distributions. We find that these quantities for the flat KPZ interfaces are reproduced by 
the 2-step RP model reasonably well, similarly to those for the circular interfaces found in agreement with the 
standard RP. Aging of the recurrence-time distributions is also clear in both cases. Note however that, while 
pdf(Tf) TjT® is known to hold for the standard RP [Eq. (|4]l] with 1 < 0 < 2 ifT^ . our 2-step RP rather indi- 

cates pdf(Tf) ^ T^ * [dotted lines in Pig.|6l[a,b)]. Since pdf(Tf) is given by the derivative of the persistence 

probability, this implies 0^ = 0±’\ hence asymptotically 0^ = 1.35 and 0'_ — 1.85 are expected for the flat 
KPZ subclass. 

In contrast to this agreement in the recurrence-time distributions, the distribution of the time-averaged sign 
a = (l/T) Jq o{x,t)dt turns out to be different between the flat KPZ subclass and the 2-step RP [Pig. Eta)], 
analogously to the results for the circular interfaces. More specifically, both the flat KPZ subclass and the 
2-step RP are found to show asymptotic broad distributions [Pig. Eta)], hence both of them exhibit WEB, but 
the distributions are again clearly different between the two systems. Note here that the time-averaged sign 
distribution for the standard RP [Eq. (|4ll] with 0 > 1 becomes infinitely narrow in the limit t ^ lITSl : this 
is however not the case here, despite 0^ > 1. The existence of the broad distribution results from the aging of 
the waiting-time distribution, i.e., from the fact that the crossover time in the waiting-time distribution grows 
with f [see Pig. |5] and Eq. (ITOl. 

The difference between the flat KPZ subclass and the 2-step RP is also detected in the correlation func¬ 
tion of sign, Csign(t,fo) = (ct(-T,f)ct(-TTo))- while our simulations of the 2-step RP show Csign(t,fo) ^ 

[Fig- Etb)], for the flat interfaces it decays as with A = 1 [Pig. Etc,d)], the characteristic exponent for 
the decorrelation of the flat KPZ subclass [see Eq. (El]. Similarly to the circular case, this difference re¬ 
sults from correlation of waiting times, which can be characterized by the generalized persistence probability 
P±,n{to + At,to). Por the 2-step RP, i.e., in the absence of correlation, we numerically find [Pig. |8|a,b)] 


P±,n>i{to + Atjo) 


^f-2e+i fQj. ^ 

At~^± for At ^ to. 


( 11 ) 


where in the latter case the two double signs are set to be the same sign for even n and the opposite ones for odd 
n. This long-time behavior can also be seen in the flat KPZ subclass [Pig.[8|c,d) for the EC experiment and (e,f) 
for the dPNG model]. In contrast, short-time behavior of P±,„>i is found to be different between the 2-step 
RP and the flat KPZ subclass [compare data and the dashed lines in Pig. [8|c-f)], the latter carrying heavier 
weight in the short-time regime. Analogously to the circular case, such pronounced short-time behavior of 
P±,n>\ seems to generate, via Eq. (El, the characteristic decay of the correlation function Csign(f,to) slower 
than that of the 2-step RP [Pig. Etb)]. 


® Strictly, since the evolution of the 2-step RP model [Eq. Gil] depends on f, it is not in the scope of the models considered 
in the renewal theory. 
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Fig. 7 Time-averaged sign distribution (a) and sign correlation function (b-d) for the flat interfaces (a,c,d) and the 2-step RP 
(b). (a) Pdf of the time-averaged sign ct = (1/T) fg cj(x,t)dr for the LC experiment (T = 15s, 30s,63s for circles, squares, 
and diamonds, respectively) and the dPNG model (T = 10000; turquoise line), compared to numerical data for the 2-step RP 
(T = 10^; dashed line). The gray vertical line indicates the ensemble-averaged value (cr) = —0.0316. The existence of the broad 
asymptotic distribution is a direct evidence of WEB in the flat KPZ subclass, but in the form different from that of the circular 
case, (b-d) Correlation function of sign, Csign(t,?o) = (o’(x,t)cj(x,?o)), at different tg for the 2-step RP (Ttot = 10®) (b), the LC 
flat interfaces (c), and the dPNC model (d). The dashed lines in the panels (c,d) indicate the exponent —A = —1, while the dotted 
line in the panel (b) shows — 6^ = —1.2. 


4 Concluding remarks 

We have shown an unexpected similarity between sign renewals of the KPZ-class fluctuations and RP, studied 
in the context of aging phenomena and WEB. Despite the fundamental difference between the two systems, 
we found, for the circular interfaces, that the KPZ waiting times obey simple power-law distributions identical 
to those defining RP, while those for the flat interfaces correspond to its straightforward extension with two 
power-law regimes [Eq. (llOb . the 2-step RP model]. Eurther quantitative agreement has been found in the 
recurrence-time distributions (Eigs.|2]and|6ll, from which the agreement in the persistence probability follows. 
These quantities have remained theoretically intractable for KPZ, but now, following the agreement we found, 
their precise forms are revealed for the circular interfaces, thanks to the exact solutions for the original RP. 
This also implies that recurrence-time statistics may be determined independently of the intercorrelation of 
waiting times, contrary to the usual beliefs. 

The correlated waiting times of KPZ otherwise generate characteristic aging properties of the KPZ-class 
fluctuations (Pigs. |3] and |7]l, especially their broad asymptotic distributions of the time-averaged sign. This 
indicates WEB of the KPZ-class fluctuations, which turned out to be different from that of RP, and in fact 
also from other types of WEB, known from the studies of single-particle observations. We therefore consider 
that the WEB found in this study is of a new kind, characteristic of many-body problems governed by the 
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Fig. 8 Generalized persistence probabilities P±_„(f,?o) for the 2-step RP (Ptot = 10^) (a,b) and for the flat KPZ-class interfaces 
[LC experiment (c,d) and dPNG model (e,f)], measured for the positive (a,c,e) and negative (b,d,f) signs. The dashed lines in all 
panels indicate the exponent —20 -I-1 = —0.6 found in the short-time regime (At ^ to) of the 2-step RP. The dotted and dot- 
dashed lines are guides for the eyes indicating exponents —0^ and —9'_, respectively, which characterize the long-time regime 
for the 2-step RP [see Eq. CD]. The same set of colors/symbols and tt is used in all panels. 


KPZ universality class. This also implies that RP cannot be a proxy for the full KPZ dymamics; instead RP 
reproduces only some of the time-correlation properties of KPZ, surprisingly well, though. 

In fact, such a partial similarity to RP was also argued in the past for the fractional Brownian motion 
(FBM), in the context of linear growth processes. Krug et al. Il28]l showed, for the stationary state of linear 
growth processes, that the stochastic process h(x,t) —h{x,tQ) is equivalent to FBM. Its first-return time (cor¬ 
responding to the waiting time of its sign) is then characterized by a power-law distribution with exponent 
6 — I— P II19II10II28II with j3 being the growth exponent, or the Flurst exponent of FBM. Cakir et al. 0 then 
suggested that the sign of FBM would form RP, showing numerical observations of its persistence probability 
as a partial support, but it turned out later that the two models behave differently in other statistical quantities 
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ini, because of the intercorrelation of waiting times. In our contribution, we studied the growth regime of 
nonlinear growth processes in the KPZ class and compared the sign of the stochastic process h{x, t) — {h{x, t)) 
with RP. As already summarized, we showed thereby precise agreement in the waiting-time distribution and 
the persistence probability, but not in the other statistical properties we studied. Understanding the mecha¬ 
nism of this partial agreement is an important issue left for future studies, all the more because no theoretical 
understanding has been made so far on persistence properties of the KPZ growth regime 12 . Such develop¬ 
ments will also help to understand the deviations from RP, which we believe carry characteristic information 
of underlying growth processes (recall our results on the correlation function). We hope this direction of 
analysis may afford a clue to elucidate hitherto unexplained time-correlation properties of the KPZ class. We 
also believe that our approach may be useful to characterize other scale-invariant processes such as critical 
phenomena. 
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Appendix A Studied systems 

In this appendix we briefly describe the three systems studied in this paper, namely the LC experiment I49I48I . the off-lattice 
Eden model (45), and t he dPN G model, all known to be in the KPZ class. The experimental results are obtained from the raw 
data acquired in Refs. I49l48l . The readers are referred to these publications for the complete description of the experimental 
system. 


Appendix A. 1 LC experiment 

The experiment concerns fluctuating interfaces between two turbulent regim es of elect rically driven nematic liquid crystal, 
called the dynamic scattering modes 1 and 2 (DSMl and DSM2, respectively) 1471491481 . The DSM1/DSM2 configuration can 
be argued to lie in pure two dimensions Q so the interfaces in between are one-dimensional. Under sufficiently high applied 
voltage, here 26 V, DSM2 is more stable than DSMl, and the interfaces grow until the whole system is occupied by DSM2. 
The initial nucleus of the DSM2 state can be introduced by shooting laser pulses. This allows us to study both circular and flat 
growing interfaces: circular interfaces grow from a point nucleus generated by focused laser pulses, while flat interfaces originate 
from a linear region of DSM2, created by linearly expanded laser pulses. 

In Refs. I49I48I . Takeuchi and Sano measured 955 circular interfaces over time length 30.5 s and 1128 flat interfaces over 
63 s, and found the characteristic statistical properties of the circular and flat KPZ subclasses, respectively. In the present study, 
we employ the same data sets which are guaranteed to belong to these subclasses, and analyze the sign of the height fluctua¬ 
tions as explained in the main text. The sign renewals are detected at every 0.5 s and 0.35 s for the circular and flat interfaces, 
respectively. 


Appendix A.2 Off-lattice Eden model 

Numerical data for circular interfaces are obtained with the off-lattice Eden model, the version introduced in Ref. (45) which is 
sometimes called the off-lattice Eden D model. While detailed descriptions can be found in Ref (45l, in this model, one starts 
with a round particle of unit diameter placed at the origin of two-dimensional continuous space. At each time step, one randomly 
chooses one of the N existing particles, and attempts to put an identical particle next to it in a direction randomly chosen from 
the range [0,2;r). If the new particle does not overlap any existing particles, it is added as attempted, otherwise the particle is 
discarded. Time t is then increased by l/N, whether the attempt is adopted or not. Particles without enough adjacent space, to 
which no particle can be added any more, are labelled inactive and excluded from the particle counter N (but can still block 
new particles). Since we are interested in the interface, or specifically the outermost closed loop of adjacent particles, particles 
surrounded by the interface are also marked inactive and treated likewise. This model was previously shown to belong to the 
circular KPZ subclass (45l The data presented in the present paper are newly obtained from 5000 independent simulations of 
time length 5000, and the sign renewals are detected at every time unit. 


^ This is because the DSM2 state, known to consist of densely entangled topological defects, needs to break surface anchoring 
on the bottom and top plates to be sustained in the system. 
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Appendix A.3 dPNG model 

Numerical simulations of flat interfaces are performed with the dPNG model. This is a discretized version of the PNG model, 
which is one of the exactly solvable models in the (1 + 1)-dimensional KPZ class OH- The evolution of the height variable 
h{x = iAx, t = nAt) of the dPNG model, with non-negative integers h, i, n, is given by the following equation: 

h{x,t + At) = mdLx{h{x — Ax,t),h{x,t),h{x + Ax,t)} + rj{x,t), (12) 

where r\{x,t) is an independent and identically distributed random variable generated from the geometric distribution, Prob(T) = 
^:) = (1 — p)'‘p with p = pAxAt. The original PNG model with nucleation rate p and nucleus expansion rate Ax/At is retrieved 
by the continuum limit Ax^O and At 0. 

In our study, we set p = 2, Ax = At = 0.1 and the periodic boundary condition h(L, t) = h{0, t) with L=l(f (or 10^ lattice 
units). We start from the flat initial condition h(x,0) = 0 and evolve the system until t = lO"* hy 10^ independent simulations. The 
sign renewals are detected at every time step, i.e.. At = Q.\ time unit. Note that the dPNG model with Ax = At = 0.\ shows the 
same universal statistical properties as the original PNG model, provided that the height variable h(x,t) is appropriately rescaled 
using non-universal scaling coefficients [Voo and F in Eq. Q]. The values of the scaling coefficients depend on Ax and At: for 
example, they are estimated at fs; 2.2408 and F fs; 1.571 for the dPNG model studied here (the evaluation method described 
in Ref. 1481 is used), while the values for the original PNG model (corresponding to Ax, At 0) are v.o = 2 and F = 1. 


Appendix B Generalized persistence probability for RP 


Here we derive two asymptotic behaviors of the generalized persistence probability for the renewal process with a power-law 
waiting-time distribution. We assume Eq. (|4j with 6 < 1 for the waiting time distribution. Thus, the Laplace transform of the 
probability density function (pdf) of waiting times t, p(t) = p'(t), is given by 

p(s)=\-as^+ e{s), (13) 

witha = r(l-0)T® QH. The generalized persistent probability can be represented by 

P„{toFAt,to)=VrQh[At„<At\to\-¥roh[Atn+\ <4f;to], (14) 

where Prob[zit„ < Afpo] is the prohahility that At„ = Tf(to) + ^2 H-h < At holds, with waiting times t, and the forward 

recurrence time 'rf(fo) (time elapsed from to to the first renewal event since then). Eor tt > 1, the double Laplace transform of 
P„(to + At,to) with respect to to and At can be calculated as follows: 


P„(s,u)= f f P„{to + At,to)e “'^'dtodAt 
Jo Jo 

u u ’ 

where fE(u;s) is the double Laplace transform of pdf(Tf(fo);to). given by ITSl 

fE{u-s)=[ [ pdf('i:f;fo)e“"^f“’'"d'rfdfo 

Jo Jo 

^ 1 

s — u 1 — p(s) 


(15) 


(16) 


Therefore, 


p{u)-p{s) 


1 




I-P(w) 


’ s — u I — p{s)' 

Now we consider the following two asymptotic limits. Eor u s (xo ^ to <t. At), we obtain 

P„{s,u) ~ 

s 


(17) 


(18) 


where we used approximation ^(u)" * ~ 1 —a(n — 1)m® ~ 1. In other words, u is so small that a{n — l)i<® ^ 1, i.e.. At 2> 
(n — l)*/®To — Then the inverse Laplace transform yields 


P„{to + At,to) 



(19) 


for Xo<^to <t. At and At ^ n^/^Xg. We note that this asymptotic behavior does not depend on tg nor n. In contrast, for ^ ^ m ^ 
T(7' (To <fo), 

Pn(s,u) ~ 


( 20 ) 
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where we used the same approximation as in the previous case. The inverse Laplace transform then yields 


P„(to + — 


r(i-0) /to 

r(0)r(2-2e) Uo 




\'^o/ 


for to- This asymptotic behavior is independent of n but does depends on to- 


( 21 ) 
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